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^ \ We discuss the gauge parameter dependence of particle spectra in statistical 

Oh" quantum electrodynamics and conclude that the electron spectrum is gauge- 

parameter dependent. The physical spectrum being obtained in the Landau 
CN ■ gauge, which leads to gauge invariance in a restricted class of gauge transfor- 

mations. We compute the thermal self-energy of electrons in a dense media 
to first order in the fine structure constant. In the zero-temperature limit and 
at high density the dispersion equation is solved for some gauges, the curves 
^NJ ■ differing from that corresponding to the physical transverse gauge, 
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I. INTRODUCTION 



The investigation of relativistic quantum electrodynamics at finite temperature (T) and 
chemical potential (/i) is of relevance in several astrophysical and cosmological contexts [|T|. 
Interest in finite temperature and density of quantum field theories has increased contin- 
uosly since the early seventies ||. The problem about the dependence of the fermionic 
Green functions on the gauge parameter is not new. For QED, it was solved by Landau and 
Khalatnikov in H and simultaneously by Fradkin in 0]. These authors showed that the 
longitudinal part of gauge field appears as a phase factor into the electron Green function 
(see also ||). The last result was extended by Fradkin || to the statistical quantum elec- 
trodynamics case, where the propagators also depend on the four-velocity of the medium. 
It is a well established fact that in QFT the electron propagator is gauge invariant on the 
tree-level mass shell (P 2 = m 2 ), see [BJ. For finite temperature, this invariance has been also 
formally shown |J. In the literature it has been explicitly verified this conclusion in several 
limiting cases |TU|-[T^]. However, it should be noted that the calculation of the mass operator 
is carried out precisely to investigate the system perturbatively out of the mass shell , where 
the condition of gauge parameter independence frequently mentioned (see i.e. ||) is not 
fulfilled, since the dispersion equation differs in general from the form P 2 = m 2 , although it 
is expected that the gauge independence of the dispersion equation remains valid near the 
mass shell. But the departure of the spectrum from the P 2 = m 2 form in the temperature 
case [^3j excludes the validity of the above mentioned proof at any order of perturbation 



theory. We will return to this point along the paper. 

A common practice in the literature has been to admit gauge independence, to choose the 
Feynman gauge and to check it perturbatively in some limits. Despite of the simplifications 
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introduced by the Feynman gauge, the Landau gauge has a special importance because the 
longitudinal part of the photon propagator contributes as a phase factor to the electron 
Green function (see for a general discussion) and thus it does not contribute to the 
physical spectrum when the gauge parameter is chosen in a such a way that the longitudial 
part of the photon propagator is zero. The purpose of the present paper is to discuss the 
problem of the gauge-parameter dependence in a non-perturbative finite temperature field 
theory, to investigate the behavior of the electron dispersion equation at finite temperature 
and density and to consider the limit of T = and high density afterwards. 

This paper is organized as follows: in section 2 we discuss the problem of the gauge- 
parameter dependence in a non-perturbative finite temperature field theory. In section 3, we 
compute the electron self-energy in the one loop approximation. In section 4 two dispersion 
equations are obtained in the limit of zero temperature for ultrahigh density. Section 5 
contain the discussion and conclusions. In appendix A we derive the scalar quantities in 
terms of which the mass operator is written, we give the definitions of the main integrals that 
appear in the calculation, and we carry out the Matsubara sums in the fourth momentum 
component k^. 



In what follows we want to argue in a perturbative-independent way that the physical 
electron spectrum can be obtained (in a covariant way) only by using the so-called Landau 
gauge 1/a = 0, which is actually compatible with a large class of gauge transformations, 
as indicated below. Usually the gauge parameter dependence a is introduced in QED to 
make it possible to invert the (otherwise singular) inverse photon propagator which is 
four- dimensional transverse. But in QFT as well as in statistics, the photon spectrum is 
a-independent. In the case of a medium (/i ^ 0, T ^ ) due to the gauge invariance of 
the theory || we can write the inverse photon Green function, by introducing in the QED 
Lagrangian the gauge ad fl A fl d u A u , as 



II. GAUGE-PARAMETER DEPENDENCE 
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and Ii^ v has the general four-dimensional tensor structure 
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The photon Green function is then 
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where 
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1 B 

Al = Bl = (k*-A)(k*-A + B(l-jgp)y W 

In the rest frame the polarization operator IT^ has the tensor structure: 

f (<% - ^) A ( k > h) + n 44 ^|| i,i = 1,2,3 
= ^ (5) 

where Il 44 = j£ ( A + ). We have also in the rest frame, in the case a — > oo, 

f fe-^)C+^§D 44 z,j = 1,2,3 

^ = < ^ (6) 

[ -D 4 i = Di4 = ^-044, 

where C = l/(k 2 - A), D M = k 4 /[k 4 {k 2 - II 44 )]. 

The photon spectrum may be obtained either by solving (|l|), in which a^O factorizes, or 
from the poles of the four-dimensional transverse part of D^. In the first case the following 
non-perturbative gauge-parameter independent and gauge invariant spectrum is obtained in 
the rest frame 

k 2 + A(k 2 ,h) = } 1+ n 44(k 2 ,A; 4 ) =o (7) 

The first equation corresponds to the spatial transverse modes and the second (after multi- 
plication by k 2 ^ 0), to the spatial longitudinal one. These are the physical modes, (which 
in the one- loop approximation and at high temperature were studied in detail by |H| ). A 
factor A; 4 , accounting for the unphysical modes introduced by the gauge condition, when cal- 
culating [ DetD^ ]~ 1//2 , is removed by the Faddeev- Popov ghost term. The photon spectrum 
can be also obtained directly from the poles of the Green function (01). We have also the 
non-physical modes k 4 = 0, given by the longitudinal term, which are not present in (P) 
since they were removed by choosing the gauge parameter a —>■ oo. In calculating the exact 
polarization operator in quantum field theory as well as in statistics, only the transverse 
electron Green function contributes. The property stems from the gauge invariance and 
gauge parameter independence of H^ u , as can be derived from the Ward identities obtained 
from the QED effective action T in vacuum as well as in the temperature case, 

^M — = ad%S( Xl - x 2 ) (8) 



which in momentum space reads 



k„D-l(k) = ak u k 2 , (9) 



where the substitution of ([[]) in (||) leads to the four- dimensional transversality for IL^. This 
means that Il^{x,y) is gauge-invariant and gauge-parameter independent 0,0. Explicitly 
we have that H til/ (x, y) is given by 
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n^(x, y) = e 2 Tr J 7 M G(x, z)T u (z, y, y')G(y', x)d 4 xd A y, (10) 

where the integration in the fourth coordinate must be understood as in the interval —/3, (3 
in the temperature case, and T(z,y,y') = ^5{z — y)5{y — y') — S f-^ ■ Thus, if a gauge 
parameter dependent G(x, y\a) is chosen, cancellations must occur in the integrand in such 
a way that the integral obtained from (|To|) must be the same as that obtained by using 
G(x,yoo) = G(x,y). This means that only the physical poles of G(x,y) contribute to H^, 
and in consequence, to the transverse part of the exact photon Green function D^. 

Concerning the Fermion spectrum, the gauge parameter dependence is introduced in the 
calculation of the mass operator E in terms of G(x, y\a) and D liv {x,y\a). Now, for the 
statistical shown by Fradkin |J, the gauge parameter dependence is given by a 

similar formula than in QED. We shall write it for the one-particle electron Green function 
as, 

G(x -y\a) = G (x - y) exp { e 2 [A* (x - y) - A (0)] } (11) 

where 

(in fllTD, Go is the propagator in the Landau gauge (1/a = 0) and in our representation 
coincides with G*), the poles of G(p, a) are obtained from the poles of the Fourier transform 
of the right hand term in ([TTj) which is the convolution of Gq(x — y) with the exponential 
a-dependent factor. Due to this convolution, in any order of a perturbative expansion of 
G(p), the resulting poles are a-dependent. 

In addition to our previous arguments concerning the electron Green function dependence 
on the gauge parameter, we want to mention the following one: Let us call E* and E' 
respectively the transverse and longitudinal terms of the electron mass operator, where 

S ' = ^E / IMP + k)T u (p + k, k)^D l d% (13) 

which is obviously nonzero and gauge-parameter dependent. Then, by taking D l {k) = 1/ak 2 
we observe it contains also the contribution of the unphysical photon modes k\ + k 2 = 0. 
By using Ward identities it can be shown || that 

2 7 

[i^P, + m + E*(p)]G(p) = 1 - £ / IMP + k )^ d ' k - ( 14 ) 



We observe also from this expression that the poles of G(p) are given by the zeros of the 
term in squared brackets in flTjj) plus the poles introduced by the second term containing 
the unphysical photon modes. 



Thus, from (|12|) or ( [L4j) it is seen that the poles of the electron Green function become 
independent of the unphysical photon modes if it is taken the limit a — > oo words, by 
eliminating the contribution of these modes to E. The zero temperature QED case can be 
treated by following similar arguments. We may cite especially ref. [14] in which a very 
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interesting discussion of the electron physical modes is given in the one-loop approximation 
for E by using the Coulomb gauge. 

Thus, although there are general proofs of gauge invariance for gauge boson spectrum 
29fl , for fermions there are not (we consider there cannot be) explicit proofs of gauge pa- 



rameter independence of the one-particle spectrum given by the poles of the Green function, 
except in very specific the infrared limit considered by Abrikosov However, 

at least in the abelian case, if all calculations are made in the Landau gauge a — > oo, which 
corresponds to choosing the longitudinal part of the photon Green function as equal to zero, 
the results lead to the physical fermion spectrum, and are gauge invariant under the class 
of gauge transformations A^(x)' = A^(x) + d fJ> r](x), ip(x)' = e ev ^-ip(x) in which r)(x) = 
—d^Aix)^/ d 2 This condition excludes the unphysical photon modes from E(x, y). 

For non- abelian theories, an expression analogous to ( |H| ) has not been obtained yet. 
However, we consider the previous arguments might be useful in considering the problem of 
gauge parameter independence of the fermion spectrum in the non-abelian case. 

III. ELECTRON SELF ENERGY 

In this section we will perform the computation of the mass operator in the one loop 
approximation for a general gauge. We shall redefine the gauge parameter as £ = 1/a 
for simplicity in the forthcoming expressions. We will keep the dependence on the gauge 
parameter to allow a further comparison of the results for differents gauges. We start from 
the expression for the electron self-energy in the one loop approximation in the temperature 
formalism (we use the shifted momenta p* = Pp — i[i e 84,\) |25| 

= T^E / d 3 k^G (p* + k)TlD%{k) (15) 



e 2 



k.4 



l <W + (£-1) KKjk 2 

KlpPu + k) 2 + m 1^ V k 2 



(16) 



Due to the medium we have the additional four velocity vector Up (in the rest frame 
m m = (0, 0, 0, i), so that k ■ u = ik± = —uj . Here, £ is a parameter of arbitrary gauge) and 
thus the vectors, on which the physical quantities may depend, are fc M and The matrix 
structure of the mass operator is now 

E (p) = i (a'jpPp + b'jpUp + c) (17) 

After some work (see the appendix), we have, for an arbitrary £ 

a = g 2 ^{(p 2 I 1 + 4I 3 + 2h + 2iuh) +£ (fh - I 3 - 2 h- 2iu I 5 )} (18) 

b = g 2 — [ [-3 i p 2 I 2 - 3 uo h + 2 uj I A - 2 i (p 2 - cj 2 ) I 5 + 2 p 2 uo 7 6 

+ £ [-ip 2 I 2 + uj J 3 + 2uj h - 2% (p 2 - uj 2 ) h + 2 p 2 to I 6 ]} (19) 
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c = (3 + £) ig 2 m e Ii 



(20) 



where g 2 = e 2 / {2tt) 2 is the fine structure constant and I\ — J 6 are integrals defined in the 
appendix. 

In the appendix, we evaluate these integrals, that is we perform the Matsubara sum and 
do some algebra, subject to the following approximations: (1) we drop the terms without 
dependence on the temperature, i.e. we ignore non-thermal electrons; (2) we assume that 
the chemical potential is much larger than the temperature, \i >> T, so that we may 
approximate the distribution function of the electrons by a step function Q (3) in accordance 
with the previous point, since the Fermi surface for positrons is negative, we assume that 
the distribution functions for positrons are exactly zero; and (4) in agreement again with all 
of this, we effectively neglect the temperature so that also the photon distribution function 
is set to zero. 

The upshot is a system with a Fermi liquid of electrons and no real photons or positrons. 
Of course, in a physical situation charge must be balanced, which in the framework we 
consider can be easily achieved through a proton background (also, the proton and electron 
tadpoles cancel each other, so we ignore tadpoles). 

From the expressions for the integrals Jj in the appendix, we carry out the analytic 
prolongation p^+i ji = iu>, the angular integrations, and finally neglect the mass J fi 2 — m 2 ~ 



/i in the upper limit in the integrals. Note that we do not neglect the mass with respect to 
the momentum. The expressions we obtain read as follows: 



h 



dk k 



In 



k 2 + 



mi 



A + B 



A-B 



(21) 



h = — 



7T I 



fdkk (uj - Jk 2 + m 2 ) A + B 
In 



V 



k 2 + m 2 



A-B 



(22) 



-2vr 



dk k 2 



k 2 + m 2 



A , 

-2 + -In 



A + B 



A-B 



P 

2~B 



In 



A + B 



A-B 



(23) 



-2vr 



dk k 2 



2p 2 (k 2 +p 2 + A) p 2 



k 2 + m 2 p [ A 2 — B 2 



2B 



In 



A + B 



A-B 



(24) 



-2%i 



dkk : 



(u — \Jk 2 + mfj 



+ ml 



A-2p 2 1 
A 2 -B 2 ~ 2B 



In 



A + B 
A-B 



(25) 



1 Actually, as soon as [i is about 10T, this approximation, which is just the characterization of the 
Fermi surface, is excellent, just as in condensed matter. 
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with the notation 



h = -2 



7T 



dkk 2 (lu — Jk 2 + m 2 



k 2 + m 2 e (A 2 - B 2 ) 



(26) 



A = p 2 — lu 2 — ml + 2cuy A; 2 + m 
B = 2kp 



(27) 



IV. DISPERSION RELATIONS 

In a medium, the propagation of an electron is described by the effective Lagrangian 
term 



L m = *e(p) S e 1 * e (p) = * e (p) (i T ^ + m e - S) * e (p) 

Accordingly, the equation of motion for the electron field is 

s; 1 * e { P ) = o 

The dispersion equation is thus 

|7 M ((1 - a) p M - ib5^ A ) - i (m e - ic)\ 2 = 



(28) 



(29) 



(30) 



By the analytic continuation ik± = —u, we get the equation for the propagation of the 
electron in the medium as 



(1 — a) 2 p 2 — (m e — i c) 2 — (1 — a) oj + b 







(31) 



In order to obtain analytic representations for the integrals Ii, we neglect the electron 
mass in the square roots: this is justified because the integrands are monotonously increasing 
in k and the upper limit of integration (~ jj) is large. At ultrahigh densities, the electron 
behaves like a massless fermionic quasiparticle. Let us illustrate the point by sketching the 
evaluation of I\ (pTT) : 



h 



71 



/ dk In 

p 



p 2 - J 2 -ml + 2k (lu +p) 



p 2 — uj 2 — ml + 2k (uj — p) 



(32) 



Introduce the short-hand 



A — 2 2 2 

A = p — uj — m e 



whereby 



h 1 

/! = --/ dk [ln|A + 2A; (w+p)| -ln|A + 2A; (u - p)\] 
V i, 



(33) 
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Integrate by parts to obtain 

77 A 

h 



P 



ln|A + 2// (u+p)\ 



A 



+yU In 



2 (a; + p) 

A + 2 li (to + p) 



2 (w-p) 



ln|A + 2/i (w -p)| 



A + 2 // (a; - p) 



— p z 



(34) 



Now, take into account only the highest orders in /i (that is, assume p small, p « p) 
to obtain finally 



Similarly, we find 



i x fa —ir— In 
P 



+ p 



10 — p 



+ 



(35) 



T -V 1 

7 2 ~ — m 

P 



uj +p 



uj — p 



+ 



(36) 



7 3 « 7r p /i In 



+ p 



u; — p 



+ 



(37) 



2 cj 2 



(38) 



4 

The analysis for J 5 « is simpler because for p —>■ we have from 



(39) 



A + 2cj£;-2p 2 



(A + 2u;fc) z -{2kp) 2 4kp 



In 



A + 2k (u +p) 



A + 2 A; (w - p) 



(40) 



1 



1 



4kp 



A + 2uk Akp \A + 2uk 







(41) 



we obtain, for the ultrahigh density limit ( /i » m e ), the following expresions for the 
functions a, b y c: 



a = g tc 



(l-0^-(3 + 20^1n 



V 



uj +p 



uj — p 



(42) 



g 2 ir 



(3 + «£l» 



w + p 



— p 



^-^ln 

2a; p 



+ p 



— p 



(43) 
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c = — (3 + £) g 2 irim e — In 

p 



u; + p 



to — p 



(44) 



Replacing these expressions for the integrals in the dispersion equation (|31|), is easy to 
obtain the dispersion relations, which constitute our main result: 



P 



> 1 9,2 A* 

z + TOg + 0. 7T — 



2p 



;i-0 — + (3 + 201n 



cj + + p 



w + - p 



(45) 



Here, uo + is the normal dispersion relation and corresponds to an electron-like quasi- 
particle. There is also a abnormal dispersion branch (c<j_ < 0), which have an essential 
collective character: 



p2 _j_ m < 



1 2 A 4 

2p 



1-0 — + (3 + 201n 



w_ +p 



— p 



(46) 



The solution of (|45|) and dig) are showed in Fig.|T] together with the light cone and the free 
dispersion law for two differents values of density. As can be observed the behavior of the 
dispersion curves is analogous to the reported in ||16|| , although in our case lower densities 
are employed, which allows us to have a more realistic description of the astrophysical and 
cosmological contexs. In Figs.[|and |3] are plotted the normal and abnormal dispersion curves 
for three values of the density. In the case of u + it is seen that its effective mass grows and 
its approach to the light cone is faster with increasing density. In the case of u>- ,its curve 
tends to approach the light cone faster than uj + , moreover it can be observed the occurence 
of a minimum for a finite value of p which decreases with the density. The abnormal branch 
shows a negative effective mass near the p = as the hole quasiparticles. 

Note that we get two different dispersion relations for particles and holes, unlike other 
results at finite temperature, such as densities near the electron mass |l7j , high temperature 
but without chemical potential [[It]], and others ]Tj|. By taking into account the electron 
mass, we correct the pathological behavior of the derivative of the dispersion curve at p = 



obtained by several authors for massless fermions [ITJr9l|I]jr2,pD . 



V. DISCUSSION 

We return to the problem of gauge invariance. As shown in Figs||, ||, the spectrum is 
dependent on the gauge parameter, but for p/m e > 1 and near the mass shell, the dispersion 
curves obtained, for some values of the gauge parameter, approach among themselves and 
to the mass shell and become gauge- independent. This in some way verifies the mass shell 
P 2 = m 2 gauge independence pointed out by several authors, but we observe it occurs for 
large momenta, i.e. in the light cone region. But for some other gauges the behavior is 
quite different. In fl5|) one can see that in the case of the gauge £ = — 3 the dispersion 
curve behaves in a drastic different way as the cases £ = 1 and £ = in the same region of 
momenta. One must stress at this point that the gauge-dependent scalar b plays an essential 
role in producing a departure of the spectrum from the P 2 = m 2 mass shell, which manifests 
especially in the low momentum limit: lim p ^ b = b(u, p, £) 7^ 0. Thus, the mass shell gauge 
invariance of [0 d° es n °t ^ ea< ^ ^° a gauge-independent spectrum. From the above results 
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we conclude that the information extracted from the curves in the region far from the mass 
shell, e.g. the effective mass, are not physical but a gauge artifact. However the results of 
section 2 indicate that the Landau gauge, in despite of introducing algebraic complications, 
is the appropriate. 

The minimum for po ^ in the abnormal solution, suggests an analogy with a superffuid 
behavior, (or superconductivity, since we are considering charged fermionics particles), but 
this conclusion can't obtained only of mass operator. To give a criterium of superconduc- 
tivity is necessary to analyze the vertex part. The electrons at high density constitute a 
weakly interacting attractive Fermi system. Obviously, to reach the conclusion that this 
system exhibits superfluidity, it is better to shift the energy axis to the point m e + 5m(/x), 
where 5m(p) = lim (uo — m e ). There, the positron-like quasiparticle spectrum is tangent at 



p = (as in the non-relativistic case |2Tfl). 

On the other hand, the tangent to the positron-like quasiparticle curve gives the velocity 
below which a superffuid effect can be expected. 

This astonishing behavior of the dispersion curve for the positron-like quasiparticle, in 
a weakly coupled Fermi gas, was already observed in a very different context by Weldon 
in 1989 [TjJ, who considered a quark-gluon plasma at very high temperature and without 
chemical potential. There, superfluidity is not intuitive at all. We think that the analogy 
with cold helium can be established with more sense in our case, of degenerate fermion gas 
at low temperature, since the phenomena of superconductivity and superfluidity disappear 
at high temperatures. 
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VII. APPENDIX 

To obtain the constants on which the mass operator depends, we proceed as in reference 
Q. First, define the auxiliar functions 

A = \ Tr [E (P) • 7%] (47) 

B = \ Tr [E (P) ■ 7%] (48) 

C = \ Tr [E (P)] (49) 

Using flT7|), we find 

a = A + iplB (50) 
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b = ±(iplA-p* 2 B) (51) 

c = C (52) 
The main integrals that appear in the calculation of these constants are 

h = VJ (k 2 + k 2 )(k 2 + k 2 + 2kv + 2k A p A + m 2 + v 2 + v 2 ) d k (53) 



f k A 

h = ^ ' (k 2 + kf) (k 2 + k 2 + 2kp + 2 k A p A + m 2 + p 2 + p 2 )^ (54) 



k,4 



J _ /* /rr\ 

3 ^ I (k 2 + k 2 A ){k 2 + k 2 A + 2kp + 2k A p A + m 2 +p 2 +p 2 A ) 1 ' 



&4 



h = E / 2 — ^ (56) 

V •/ (k 2 + k 2 f (k 2 + k 2 + 2 kp + 2 k A p A + m 2 + p 2 + p 2 ) 

h = E / 2 — d 3 k (57) 

V J {k 2 + k 2 A f {k 2 + k 2 + 2 fc p + 2 k A p A + m 2 + p 2 + p 2 ) 

h = E / 2 ^ (58) 

{k 2 + k 2 f{k 2 + k 2 + 2kp + 2k A p A + m 2 +p 2 +p 2 A ) 

The integration is over A;-three-space and the Matsubara sum runs over the discrete fourth 
component of momentum, which for bosons is taken as k A = 2niij3 (n 6 Z). 

Assuming /j, » T, only thermal electrons survive with a step-function distribution, 
whereas both positrons and photons get killed, and we find 

h = I 7T. f ~< 9 ~ -< ^ k (59) 

J 2ie k+p [-p A -i{e k - e k+p + [-p A + i (e k + e k+p - fj,)\ 

I 2 = [ = [-P* + *^-^H»-e k+P ) 3k 

J 2ie k+p [-p 4 -i{e k - e k+p + n)\ \-p A + i{e k + e k+p - n)\ 



f k -l?0(n-£k +P ) d 3 k (61) 

J 2ie k+p [-p A -i(e k - e k+p + //)] \-p A + i{e k + e k+p - //)] 



f (Z-?)e(»- e k+p ) ^ 

J 2ie k+p [(p A + iji) +i(e k - e k+p )f[(p A + iji) -i(e k + e k+p )} 2 
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h = - 



k [-(p4 + ifi)+i e k +p] 0(y- gfc+p) 

2ie k+p [(p 4 + ifi)+i (e k - e k+p )] 2 [(p 4 + i //) - % (e k + e k+p )f 



(63) 



/ 2% 



[- (p4 + in)+i e k+p f 9(fi- e k+p ) 



2ie k+p [(p 4 + i //) + i (e k - e k+p )f [(p 4 + i //) - i (e fc + £fc+ P )] 2 



cfA; 



(64) 



In the above, e k+p = \J {k + p) 2 + m 2 , 



fc 

function characteristic of the low temperature limit. 



and 6 (A) is the Heaviside unitary step 
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FIGURES 




FIG. 1. The spectrum of fermionic excitations in an ultrarelativistic plasma for two values of 
density. In dots are showed the light cone and the free dispersion relation. 
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FIG. 2. Dependence of the normal branch with density. The dots lines correspond to the light 
cone and the free dispersion relation. 
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FIG. 3. Dependence of the abnormal branch with density. The dots lines correspond to the 
light cone and the free dispersion relation. 
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FIG. 4. Dependence of the normal branch with the gauge parameter. The dots and dashed 
lines correspond to the light cone and the free dispersion relation respectively. 
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